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The quark generalised parton distributions of the pion are calculated from lightcone 
wavef unctions in transverse lattice gauge theory at large N^.. The pion effective size is found 
to decrease with increasing momentum transfer. An analytic ansatz, consistent with finite 
boundstate lightcone energy conditions, is given for the lightcone momentum dependence 
of the wavef unctions. This leads to simple, universal predictions for the behaviour of the 
distributions near the endpoints, complementing numerical DLCQ data. 
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I. INTRODUCTION 



The generalised parton distribution of the pion has been studied in a number of different 
models Following earlier work 0, van de Sande and the author Q recently computed 
numerically the lightcone wavefunctions of the pion in transverse lattice gauge theory at 
large Nc 0,1^ (see also [^). In this talk, these wavefunctions will be applied to deduce the 
generalised quark parton distributions in both transverse momentum and position space ■ 
We also present an analytic ansatz for the lightcone momentum dependence of the wavefunc- 
tions that leads to simple, universal predictions for endpoint behaviour of the distributions, 
complementing the numerical data. 

In transverse lattice gauge theory in 3 + 1 dimensions, two coordinates x", a G {0, 3}, are 
continuous while two directions x = are discrete. There are longitudinal continuum 

gauge potentials ^0,^3, transverse flux link fields Mi,M2, and fermions For transverse 
lattice spacings a of order the hadron size, the strategy is to perform a colour-dieletric 
expansion of the most general lightcone Hamiltonian, renormalisable with respect to the 
continuum coordinates, in powers of Mi, M2 after eliminating Aq, A3 by lightcone gauge 
fixing. Provided the fields \E',Mi,M2 are chosen sufficiently heavy, one can truncate this 
expansion to study the low lying hadron boundstates dominated by just a few particles of 
these fields. The remaining free couplings in the effective Hamiltonian are fixed by optimizing 
symmetries broken by the regulator and, if necessary, phenomenology. 

We take the pion wavefunctions from Ref. which contains details of the construction 
of the lightcone Hamiltonian and Fock space, the renormalisation and the determination 
of various residual couplings appearing in the lightcone hamiltonian. The couplings were 
constrained by optimizing covariance of low-lying meson and glueball wavefunctions, rota- 
tional invariance of the heavy source potential, and fitting two phenomenological parameters 
in the meson sector, conveniently taken to be /^r and the p meson mass. In addition, two 
fundamanetal scales are taken as given, the tt meson mass and the string tension ^/a. 

A typical large- A'^c meson eigenstate state consists of linear combinations of gauge invari- 
ant basis states formed from quarks separated by a chain of transverse links \I'MM ■ ■ ■ M\E'. 
If we label the transverse orientation of link fields by indices Xj G {±1,±2}, the partonic 
decomposition for a meson of momentum P = (P^, P^, P) when the transverse momentum 
P = can be written 



00 „ 

\P\P = o) = Y.J[dx]n E 



E 



1pn{Xi, h, h', Xj) X 

n=2 zi,z„,h,/i' 

|(xi,/i,zi);(x2,Ai,Z2));---;( 

where h, h' denote quark and anti-quark helicities respectively, Zj is the transverse position 
and Xi the fraction of P^ carried by the z*^ parton. We have 



j [dx\n = J dxi - ■■ dxnS I^E ^' " 1^ 



(2) 



and 



indicates that the sum over orientations of links must form an unbroken chain on 
the transverse lattice between quark and anti-quark. We define a set of hadron states boosted 
to general transverse momentum P by applying the Poincare generators M+ = (M+^, M+^). 
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This gives for each parton Fock state 



exp 



-iM+.P/P+ \{xi,h,z 



exp [iP.c] \ {xi,h,zi] 



(3) 



^ ] XiZii 

i=l 



(4) 



The hghtcone wavefunctions ipnixi, h, h', Xj) computed in Ref.fs'l were subject to a number 
of cutoffs. Firstly, the transverse lattice spacing is held fixed at a ~ 300MeV~^, representing 
the renormalisation scale. Thus, it would be inappropriate to use momentum transfers Q 
much larger than this. A cutoff n < 5 was used. This translates into a maximum transverse 
size for the meson of 3 links, or about 2 fm. The effect is to reduce slightly the charge radius 
of the pion from its experimental value, as was shown in the form factor -F(Q^) calculated 
in Ref. [7|. Finally, calculations were performed in DLCQ, which discretizes the momentum 
fractions Xi in units of a small cutoff 1/K, to create a finite-dimensional Fock space. The 
largest basis studied was K = 20, which is of dimension 82,470. Observables tend not to 
change much when extrapolated from K = 20 to = oo, except close to the edges of phase 
space. We will give an analytic approach to these edges in section HVl 



II. NUMERICAL GENERALISED PARTON DISTRIBUTION 



We follow the conventions of Diehl [9|], defining the generalised distributions in lightcone 

gauge as 



n{x,^,Q' 



1 



4n 



(Pin|^(-;2-/2)7+^(^-/2)|Po,t) 



(5) 



Q = Pi, - Pout (6) 

(-Pin ~ -Pout)^ 

2P^ 

(-Pin + -Pout)^ 



p _ (-Pin -Pput)^ /y\ 

^ " 2P+ 



(8) 

Unfortunately, we are not able to explore the region a; < ^ with the strict large-A^^c 
wavefunctions. That region is governed by their l/N^ corrections, since it depends upon 
quark pair production. Even at a; > ^, in DLCQ one needs to form overlaps of wavefunctions 
with different K in fixed ratio for fixed ^. Existing data do not cover enough values of K. 
Thus, we concentrate on ^ = for the DLCQ analysis. 

The i^' = 20 data are displayed in fig.l for two extreme (!) values of Q at ^ = 0. T-C{x, 0, 0) 
is the conventional structure function which exhibits rising Regge behaviour at small x 
characteristic of a system with gluonic degrees of freedom. It falls to zero precipitously near 
a; — > 1. This may be an indication that the chiral limit is near (the pion mass was fit to 
140 MeV), since chirally symmetric quark models 0] produce a flat distribution. As purely 
transverse Q is increased, the distribution is seen to deplete at small x but is Q-independent 
as X — i> 1. 
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FIG. 1: Generalised quark distribution of the pion at zero skewdness and DLCQ cutoff K = 20: 
black points are Q = 0, grey are Q = 500 MeV. 
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FIG. 2: Effective size of tlie pion as measured by Eq. Q. 



A measure of the effective size of the pion is given by the average distance between the 
active quark and the center of momentum of the spectators ^ij 



d 



F{Q^) dQ^ Jo 



XI 



(9) 



This is plotted in fig. 2. The authors of Ref. [uj claim that B^(Q^ oo) — > is a condition 
for colour transparency. The data are consistent with the pion size reducing, but obviously 
it is impossible to reliably extrapolate into the large momentum transfer region. 
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FIG. 3: Impact parameter dependent valence quark distributions of the pion. (a) The fuh distri- 
bution TC{x, 0,0), which sums over ah impact parameters b. (b)-(e) The sum of contributions at 
impact parameters b along lattice axes such that 6 = |b| = 0, a, 2a, 3a respectively. 



III. IMPACT PARAMETER DISTRIBUTION 



The impact parameter dependent quark distribution was analysed in Ref. 0]. It is defined 

by H 



X(x,e,b) 



d^ ^ixP'^z- ^ 



^(-z-/2,b)7+vI/(^-/2,b)|P+,bn 



where 



■uout 



e 



:b , 



-b . 



(10) 



(11) 



bo is the transverse position of the hadron, assumed spinless. X(x, 0, b) is simply the prob- 
ability of a quark carrying fraction x of the lightcone momentum P"*" when at transverse 
position b [l3[. 

To translate results to hadron impact parameter space in the transverse directions, we 
integrate P over the Brillouin zone to obtain a hadron state localised at transverse position 




oo „ 

|p+,bo = o) = [dAn E 

o _ _ 



E 



a? sin c^ir/a sin c^vr/ a 



X 



n=2 zi,Zn, h,h' 

ipnixi, h, h', \j)\{xi, h,zi); ■ ■ ■ ; h', z„)) 



(12) 
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Hence 



oo „ 

J(x,0,b) = Yl j[dx]n E 

n=2 x\,Xn,h,h' 



E 



sinc^7r/asinc^7r/a 



Cl 



C2 



— x)5(zi — b) X 



\i)n{Xi,h,h\\j)\^ . 



(13) 



Due to the lattice cutoff, this should strictly be interpreted as the probability of finding 
the quark within one lattice spacing of impact parameter b. We show results for the first 
few values of 6 = |b| sampled along lattice axes in fig. 3. In this case, the curves are the 
result of an extrapolation to = oo from data dX K = 10, 12, 15, 20. Fig. 3(a) shows the 
full distribution 7i(a;, 0, 0), which sums over all impact parameters b (one can compare this 
with the if = 20 version in Fig. 1). Fig. 3(b) is the contribution from b = (0,0), Fig. 3(c) 
is the sum of contributions from b = (a, 0), (0, a), (—a, 0), (0, —a), etc.. This clearly shows 
that the valence region x > 0.5 is dominated by quarks with impact parameter less than the 
charge radius ~ 2/3 fm. It also suggests that the rapid rise at small x is dominated by larger 
impact parameters b > 2/3 fm. This all implies a sharp fall-off of the hadron wavefunction 
in transverse space at a particular (x-dependent) radius. 



IV. ANALYTIC WAVEFUNCTION ANSATZ 

The boundstate problem can be cast into a set of integral equations for the wavefunctions 

M2^„(x„/i,/i',A,) = JlC-^Pn (14) 

where, on the right hand side, the wavefunction is convoluted with a kernel /C [3] that 
exchanges lightcone momentum, creates, and annihilates partons. Typically the kernel is 
singular when one or more momentum fractions Xi vanish in this equation. Demanding that 
the boundstate mass M remains finite, one can constrain the behaviour of the wavefunctions 



ipn in these limits [1^. A thorough analysis of the endpoint behaviour of transverse lattice 
wavefunctions will be given in Ref. jl^. Here we will give the simplest ansatz for the first 
few ipn in a large- A'^c meson that has the correct endpoint behaviour when any number of Xi 
vanish: 

Mx) = C2x"(l-x)" (15) 
^^^^'^^ = ^'[{x + y)il-xW-^^y^ 

A - r - - mix + y + z){l- xr^^-'/' 

iJ,{x,y,z) - 64 [(x + y)(i-x~y)r^y^{y + zyP ^^^^ 

The exponents a and (3 are independent of n, quark helicity, and the transverse shape, 
while the coefficients C„ are dependent on all of these factors. They are all determined by 
the solution to Eqs. (fT^ . The coefficients C„ are complex in general, although overlaps of 
wavefunctions appearing in 7i are real because the theory is PT invariant. For normalis- 
abihty one must have /3 < 1/2. The simplest forms above could be generalised for use as 
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a variational basis, by distinguishing quark and link-field Xi and/or multiplying them by 
a complete basis of analytic functions of the Xj. In principle, one can use DLCQ data to 
determine the unknown parameters but, in practice, it is difficult because convergence is 
slow in the endpoint regions. We could, however, obtain a fairly reliable fit a ~ 1/2 from 

the data of Ref. Q, implying a quark distribution function of the form \Jx{l — x) at this 
low scale. 

Even without determining the unknown parameters, one can deduce a number of inter- 
esting conclusions. Firstly, one notes that the diverging singularity of the wavefunctions 
ipn becomes stronger as more momentum fractions vanish. Intuition from the free light- 
cone kinetic energy ~ Y^ii^ilxi, for partons of mass mj, would have suggested the complete 
opposite. Rather than simply vanishing at the edges of phase space, as indicated by the 
free theory, the interactions force relations between wavefunctions with different numbers 
of partons to cancel the divergences in Eq. (fT^ . This was noted for continuum QCD for a 
single vanishing Xi in the first reference of Ref. (l^J , but for the transverse lattice we are able 
to generalise it to any number of vanishing momenta. This allows us to make a number of 
definite statements relevant to generalised quark parton distributions: 

1. each Tpn contributes const, x (1 — x)^" to ^, Q"^) as x — 1; 

2. il)n contributes (In ^)"~^, > 2, to H{x, 0, 0) as x — 0; 

3. each n > 2, contributes a finite constant to ?-^(x, ^, Q'^) as x ^ ^ from above; 

4. m{x,Q,Q'^)/dQ%2^Q ~ [l-xf as X ^ 1. 

Provided the constants fall sufficiently fast with n, continuity and item 3 above suggest 
that is finite. 

V. OUTLOOK 

The transverse lattice has so far proven to be the only systematically improvable non- 
perturbative lightcone method that preserves the gauge invariance crucial to confinement. 
The current DLCQ transverse lattice calculations for mesons can be improved by going to 
higher orders of the colour dielectric expansion and testing directly for P, T, and chiral 
symmetry. As well as the initial applications to the pion presented here, one could also 
calculate the generalised distributions of the heavier mesons. Of course, most experimental 
data derive from the nucleon, so it is important to develop computations in the baryon sector. 
These have proven awkward in the past because the numerical work does not simplify in 
the \sxge-Nc limit. However, many of the analytic arguments above will also be valid for 
baryons. 
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